Cosmological constraints on extended Galileon models 
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The extended Galileon models possess tracker solutions with de Sitter attractors along which the 
I , dark energy equation of state is constant during the matter-dominated epoch, i.e. uide = — 1 — s, 

■ where s is a positive constant. Even with this phantom equation of state there are viable parameter 

' spaces in which the ghosts and Laplacian instabihties are absent. Using the observational data of 

^S) , the supernovae type la, the cosmic microwave background (CMB), and baryon acoustic oscillations, 

. ■ we place constraints on the tracker solutions at the background level and find that the parameter s 

^ ' is constrained to be s = 0.034^q q34 (95 % CL) in the flat Universe. In order to break the degeneracy 

Q between the models we also study the evolution of cosmological density perturbations relevant to 
the large-scale structure (LSS) and the Integrated-Sachs- Wolfe (ISW) effect in CMB. We show that, 
QQ ' depending on the model parameters, the LSS and the ISW effect is either positively or negatively 

correlated. It is then possible to constrain viable parameter spaces further from the observational 
. data of the ISW-LSS cross-correlation as well as from the matter power spectrum. 

u: 

' I. INTRODUCTION 

Q , The main target of the dark energy research over the next few years or so is to distinguish between the A-Cold-Dark- 
%-[ ' Matter (ACDM) model and dynamical models with time- varying equations of state wde- From the observational data 
C/3 ^ of WMAP7 combined with baryon acoustic oscillations (BAO) [1] and the Hubble constant measurement [2], Komatsu 
, et al. [3] derived the bound wde = — 1.10±0.14 (68 % CL) for the constant equation of state. Adding the supernovae 

\ type la (SN la) data provides tighter constraints on wdE) but still the phantom equation of state (wde < —1) is 
y—^ . allowed by the joint data analysis This propertypersists for the time- varying dark energy equation of state with 
^ ' the parametrization such as iude = + 'u^a(l — a) 01 1 where a is the scale factor fF]. 

I In the framework of General Relativity (GR) it is generally difficult to construct theoretically consistent models of 
dark energy which realize wde < — 1- In quintessence with a slowly varying scalar-field potential, for example, the 
field equation of state is always larger than —1. A ghost field with a negative kinetic energy leads to wde < — 1 0) 
• ' but such a field is plagued by a catastrophic instability of the vacuum associated with the spontaneous creation of 

I ghost and photon pairs @. 
I ' In modified gravitational theories it is possible to realize wde < — 1 without having ghosts and Laplacian-type 
, instabilities (see Refs. [9]). In /(i?) gravity, where the Lagrangian / is a function of the Ricci scalar R, the dark 
• • ■ energy equation of state crosses the cosmological constant boundary (wqe = -1) [IMl 

for the viable models 

constructed to satisfy cosmological and local gravity constraints [lOl - KLa [l4| . This is also the case for the Brans-Dicke 
theory [l5| with a field potential which accommodates the chameleon mechanism [l6| to suppress the propagation of 
the fifth force [T7| . In modified gravity models of dark energy based on the chameleon mechanism (including f{R) 
gravity) , the effective potential of a scalar degree of freedom needs to be carefully designed to pass cosmological and 
local gravity constraints [l8| . 

There is another class of modified gravity models of dark energy in which a nonlinear self-interaction of a scalar 
degree of freedom can lead to the recovery of GR in a local region through the Vainshtein mechanis m |19il . The 
representative models of this class are those based on the Dvali-Gabadadze-Porrati (DGP) braneworld [20| and the 
Galileon gravity [2TI| (see Refs. [l^, for the implementation of the Vainshtein mechanism in these models). The 
nonlinear interaction of the form (9(/))^n(/), which appears from the brane-bending mode in the DGP model [22,], gives 
rise to the field equation invariant under the Galilean shift 9^0 ^ -1- &p in the flat spacetime. This was extended 
to more general field Lagrangians satisfying the Galilean symmetry in the limit of the Minkowski spacetime (2ll. [23|. 

The cosmology based on the covariant Galileon or on its modified versions has been studied by many authors 
(25I [26| . In Refs. [27i [28j the dynamics of dark energy was investigated in the presence of the full covariant Galileon 
Lagrangian. In this model the solutions with different initial conditions converge to a common trajectory (tracker). 
Along the tracker the dark energy equation of state wde changes as —7/3 (radiation era) — 2 (matter era) — 1 
(de Sitter era) [2^^ - [30| . There exists a viable parameter space in which the ghosts and Laplacian instabilities are 
absent. However, the joint analysis based on the observational data of SN la, CMB, and BAO shows that the tracker 
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solution is disfavored because of the large deviation of wde from —1 during the matter era [30|, |3l|. The solutions 
that approach the tracker only at late times are allowed from the combined data analysis [3l| . 

As an extension of the covariant Galileon model, Deffayet et al. [s^l obtained the most general Lagrangian in 
scalar-tensor theories with second-order equations of motion. In four dimensions the corresponding Lagrangian is 
of the forrnJlT]) with the four functions ([^-([5 t given below. In fact this is equivalent to the Lagrangian found by 
Horndeski [33| more than 3 decades ago |34l . ISSj . The conditions for the avoidance of ghosts and Laplacian instabilities 
were recently derived in Ref. [sg'] in the presence of two perfect fluids (non-relativistic matter and radiation). 

The covariant Gahleon corresponds to the choice K = -C2X, G3 = c^X/M^, G4 = - aX'^/M^, G5 = 

3c5X^/A/^ in Eqs. ©-(H]), where Ci's are dimensionless constants, X = — 9^09u(/)/2, Mpi is the reduced Planck mass, 
and M is a constant having the dimension of mass. Kimura and Yamamoto [30] studied the model with the functions 
K = —C2X, G3 = C3M^~'^"X" (n > 1). G4 = Mpj/2, and G5 = 0, in which case the dark energy equation of state 
during the matter era is given by wue = — 1 — s with s = l/(2n— 1) > 0. At the background level this is equivalent to 
the Dvali- Turner model '37"|, which can be consistent with the observational data for n larger than the order of 1. If 
we consider the evolution of cosmological perturbations, the LSS tends to be anti-correlated with the late-time ISW 
effect. This places the tight bound on the power n, as n > 4.2 x 10'^ (95% CL) [sll, in which case the dark energy 
equation of state is practically indistinguishable from that in the ACDM model. 

In Ref. [15 the present authors proposed more general extended Galileon models with the functions K — 

-C2M^'-^'P'^XP^, G3 = C3M^-^P'XP^, G4 = - CiM^-^P^XP^, and G5 = ScsM^^^^^^'^ XP^ , where the 

masses M^'s are fixed by the Hubble parameter at the late-time de Sitter solution with (/> = constant. For the powers 
P2 — p, P3 — p + {2q — l)/2, Pi — p + 2q, P5 = p + (6g — l)/2, where p and q are positive constants, there exists 
a tracker solution characterized by Hcj)^'' = constant. During the matter-dominated epoch one has wde = ^1 — s, 
where s = p/{2q), along the tracker. This covers the model of Kimura and Yamamoto [30'| as a specific case {p — 1, 
q = n — 1/2, C4 = 0, C5 = 0). In the presence of the nonlinear field self-interactions in G4 and G5, the degeneracy 
of the background tracker solution for given values of p and q is broken by considering the evolution of cosmological 
perturbations. Hence the ISW-LSS anti-correlation found in Refs. (30l . [ssj for 04 — cc^ ~ should not be necessarily 
present for the models with non-zero values of C4 and C5 . 

In this paper we first place constraints on the tracker solution in the extended Galileon models by using the recent 
observational data of SN la, CMB, and BAO. The bound on the value s = p/{2q) is derived from the background 
cosmic expansion history with/without the cosmic curvature K. We then study the evolution of cosmological density 
perturbations in the presence of non-relativistic matter to break the degeneracy of the tracker solution at the back- 
ground level. We will show that the LSS and the ISW effect are either positively or negatively correlated, depending 
on the parameters C4 and C5. This information should be useful to distinguish between the extended Galileon models 
with different values of C4 and C5 from observations. 



II. BACKGROUND FIELD EQUATIONS 



We start with the following Lagrangian 
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i=2 



where 



C2=K{X), (2) 

£3 = ~G3{X)U4>, (3) 

£4 = G4(X) R + Gi.x [(□<^)' - (V^V.c/.) (V^V-c/.)] , (4) 

£5 = G^iX) G^, [V^V^) - (G5,x/6)[(n0)' - m<t>) (V^V,0) (V^V^,^) + 2(V'^V„0) (V"V^0) (V^V^^)] . (5) 

K and Gi {i = 3, 4, 5) are functions in terms of the field kinetic energy X = —d^4)d^(t)/2 with Gi,x = dGi/dX . R is 
the Ricci scalar, and G^^, is the Einstein tensor. If we allow the 0-dependence for the functions K and Gi as well, 
the Lagrangian ([1]) corresponds to the most general Lagrangian in scalar-tensor theories (32l . |33| . In order to discuss 
models relevant to dark energy we also take into account the perfect fluids of non-relativistic matter and radiation 
(with the Lagrangians Cm and Cr respectively) , in which case the total 4-dimensional action is given by 

j4 



/ <fx^{C + Cm + Cr). (6) 
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In the following we focus on the extended Galileon models [36j in which K and Gi are given by 



K 



-c2m: 



4(1-P2) 



XP2 



Ga 



AfpV2-C4M4 



2-4p4 



^P4 



3C5M5 



■(1+4P5) 



(7) 



where Mpi is the reduced Planck mass, ci and (i = 2,3,4,5) are dimensionless constants, and (i = 2,3,4,5) 
are constants having the dimension of mass. In the flat Universe it was shown in Ref. [36| that tracker solutions 
characterized by the condition H^P''' = constant {q > and a dot represents a derivative with respect to cosmic time 
t) are present for 



P2^P, P3^p+i2q-l)/2, 



P4=p + 2q, 



P5=P+i6q- l)/2 . 



(8) 



The covariant Galileon ^] corresponds to p = 1 and q — 1/2, i.e. P2 — Ps — 1, P4 — P5 — 2. 

We will extend the analysis to the general Friedmann-Lemaitre-Robertson- Walker (FLRW) background with the 
cosmic curvature K: 



ds^ 



dr2 



1 - Kr- 



+ r2(d6i2 +sin2 6'd02) 



(9) 



where a{t) is the scale factor. The closed, flat, and open geometries correspond to K > 0, K ~ and K < i)^ 
respectively. For the theories given by the action ^ the dynamical equations of motion are 



■iH'^Ml^ = PDE + A 
(3i72 

Pm + iHpra = , 
Pr + 4,Hpr = , 
PK + 2HpK - . 



Pr 

2H)Ml, = -Pde - Pr/i + PK/i : 



PK 



(10) 
(11) 

(12) 
(13) 
(14) 



Here H = a/a, pK 
respectively, and 



-3-R'Afpj/a^, pm and pr are the energy densities of non-relativistic matter and radiation 



PDE = 2XK^x -K + m^XGs^x - GH^Gi + 24H^X{G4,x + XGi^xx) + 2H^X{5G5.x + 2XG5,xx) , (15) 



Pbe^K- 2X4)G3,x + 2(3iJ2 + 2H)G4 - 4:{3H^X + HX + 2HX)G4,x - 8HXXG4,xx 
-2X{2H^<j) + 2HH<j) + 3H^^)G5^x - AH'^ X^'<j)Gz,xx , 

where G4 = G4 - MlJ2 ^ -c^mI'^^^XP^ 

From Eqs. ([TO)) and (fTTj) we find that there exists a de Sitter solution characterized by ij = and 
to discuss the cosmological dynamics we introduce the dimensionless variables 



(16) 



0. In order 



l+2q 



r2 



X 
XdS 



P + 2q 
l + 2q 



Pr 



(17) 



where x = (f>/{HMpi), and the subscript "dS" represents the quantities at the de Sitter solution. We relate the 
masses Mi (i = 2,--- ,5) in Eq. ^ with H^s, as M2 ee (i^ds^^pl)^^^ = (-^ds^^'M V^''')^/^^~''^"\ M4 = 



(H~^P^ M^~^P^y/P-^P^\ and M5 = (ij2+2p5^2p5-i^i/(i+4p5)^ The existence of de Sitter solutions demands that 
the coefficients C2 and C3 are related with C4 anci C5, via 



C2=o — {3a -4(3 + 2), 



C3 



V2 



2p^ 



''dS 



p+q 



— [3{p + q)ia~/3)+p] 



where 



4(2p4 - 1) 
a = 



"-dS 



P4 



C4 , 



P^2V2p,{tf 



P5 + 1/2 



C5 , 



The density parameter of dark energy, fl^E = Pbe/{3H^M^^), can be expressed as 



(18) 



(19) 



DE 



n {n [I2(a - (3){p + q)+4p- ri{2p - l)(3a - 4^ + 2)] - 3Q'(2p + iq + 1)} + 4^(p + 3q + 1) . (20) 
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From Eq. it follows that ^de + + fir + f^if = 1, where Q,^ = ^^/(aiJ^M^j) and = pk /i^H'^M^^). 

The autononious equations for ri, r2, and fir are written in terms of ri, r2, fir, a, /3, p, q. As in the case of the 
flat Universe [s^ one can show that there is a fixed point for the differential equation of ri characterized by 

ri = l, (21) 

From the definition of ri in Eq. ^T7\ this corresponds to the tracker solution where H(j)^'^ — constant. Along the 
tracker the autonomous equations for r2 and fir are 

, ip + 2q){nr+3~3r2~nK) 

^ ^^^q ' (''^ 

^ 2q{nr-l-3r2~n^)~ipr2^^ 
pr2 + 2q 

where a prime represents the derivative with respect to iV = In a. Combining these equations, we obtain the integrated 
solution 

r2=cia^(^+^)l7i+% s = f, (24) 

Zq 

where ci is a constant. For the theoretical consistency the parameter s is positive [s^. Since r2 oc 

quantity r2 grows toward the value 1 at the de Sitter solution. Along the tracker the density parameter (j20p is given 

by 

»DE = r2 = ^ ' ~ ~ e4(^+-)^ »^+- , (25) 

where the subscripts "0" represent the values today (the scale factor ao = 1, i.e. A^'o = Inao = 0). Using the relation 
nx/^r = {i^K,o/^rfi)e'^^ as weh, Eq. ^ reads 



^, _ 1 - a- + »A-.oe'^ ^^r/^^r,o + (1 - »^,o - 17^,0 - ^j^^p) (3 + 4s) e^(^+''^^n}+ynl+' ^ 



This equation can be solved as 

1 - 17,„,o - ».,0 - ^^X,0 ^4(1+.)7V f^l + s ^ ^^^^ + f], + ^ e^^ = 1 , (27) 

which is nothing but the relation 51de + + + = 1 ■ From Eq. the dark energy density parameter evolves 
as r^DE oc 7J^2(i+s) g^j-^j^ hence 



fj /O \ l/[2(l+s)] 
/ "DE.O 



Hq \ fluE 



(28) 



Since it is not generally possible to solve Eq. (|27|) for fi^ in terms of N (apart from some specific values of s such as 
s — 1), we numerically integrate Eq. (|26p and find the expression of H/Hq by using Eqs. ([25l) and ((28|). 

Along the tracker the dark energy equation of state i^de = ^de/pde and the effective equation of state Wcs — 
-1 - 2i7/(3iJ2) are given by 

3 + s(3 + - riif ) r2(3s + 3 - sfiif) - fir .„„x 

tWDE = TT— ^ , Wcff = ^TT— ^ ■ (29 

3(1 + sr2) 3(l + sr2) 

In the early cosmological epoch (r2 ^ 1) these reduce to wdb — —1 — s(3 + fir — riA')/3 and Wcs — During 
the matter era in which {ilr, \^k\} ^ 1 it follows that wde — — 1 — s < — 1 (for s > 0) and Wcs — 0. At the de Sitter 
fixed point (r2 = 1) with 17^ = = one has wue — WcS ^ —1. In Fig. [1] we plot the evolution of i«de versus 
the redshift z ~ aa/a — 1 in the open Universe with flK.o — 0.1 for the model parameters p — 1, q — 5/2, a = 3, 
P = 1.45 (i.e. s — 0.2). The tracker is shown as a solid curve, along which wde changes as —1.267 (radiation era) — > 
— 1.2 (matter era) — — 1 (de Sitter era). The effect of the cosmic curvature fix becomes important only for the late 
Universe, which affects the luminosity distance in the SN la observations. 
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log (z+1) 



Figure 1. Evolution of wde versus the redshift z for the open Universe with Q.k,o ~ 0.1 for p — 1, q = 5/2, a — 3, l3 = 1.45. 
The solid curve corresponds to the tracker solution with the initial conditions ri — 1, r2 — 1.0 x 10""^", Qr = 0.99987, 
Qk = 4.0 X 10"^^ at logio(2 + 1) = 7.245. In the cases (A) and (B) the initial conditions are chosen to be (A) ri = 1.0 x 10 , 
ra = 1.0 X 10-2^ Qr = 0.99985, = 5.0 x 10"^^ at logio(z + 1) = 7.210, and (B) n = 3.0 x 10"®, r2 = 1.0 x 10"^", 
Qr = 0.9998, = 1.15 x 10"" at log^oi^ + 1) = 6.967, respectively. 



If the solutions start from the regime ri ^ 1, the evolution of wde is different from Eq. before they reach the 
tracker. For ri ^ 1 and r2 ^ 1, wde and Wcs are approximately given by 

"°^""2(2p + 6,-l)' -<=*^-3"- (30) 

In Fig. [1] we show the variation of wde for p = 1, g = 5/2, a = 3, /3 = 1.45 with two different initial conditions 
satisfying ri ^ 1. In both cases the density parameter ilj^ today is flK,o = 0.1. The cases (A) and (B) correspond to 
the early and late trackings, respectively. For smaller initial values of ri the tracking occurs later. As estimated by 
Eq. (|30p . wde starts from the value wue — —1/16 in the deep radiation era. If the solutions do not reach the tracker 
during the matter era (as in the case (B) in Fig. [IJ, wde temporally approaches the value —1/32. 

In Ref. (3l| it was shown that the tracker for the covariant Galileon (s = 1) is disfavored from observations, but 
the late-time tracking solution is allowed from the data. This property comes from the fact that for the late-time 
tracker the deviation of wde from —1 is not significant. For s ^ 1 even the tracker is expected to be allowed from 
observations. In such cases the solutions starting from the initial conditions with ri <^ 1 should be also compatible 
with the data (because even for s = 1 the late-time tracking solution is allowed). In the following sections we will 
focus on the tracker solution to discuss the background observational constraints and the evolution of cosmological 
perturbations. 



III. OBSERVATIONAL CONSTRAINTS ON THE EXTENDED GALILEON MODELS 



In this section we place observational constraints on the tracker solution from the background cosmic expansion 
history. We use three data sets: 1) the CMB shift parameters (WMAP7) 3; 2) the BAG (SDSS7) [H; 3) and the 
SN la (Constitution) [3^. The total chi-square Xtot for all three combined data sets will be calculated on a grid 
representing a chosen set of available parameters. We then find the minimum on this grid, and consequently find the 
1(7 and 2(7 contours. 

In order to integrate Eq. (|26p . once a set of model parameters (in this case not only s, but also f2r,0: ^m.o-, and flKfl) 
is given, we have the choice of one initial condition, that is fJ^.i = ^r{Ni). In principle it is possible to solve Eq. (|26p 
backwards for a given value of flr{0) = ^r,o, but we find that the integrated results are prone to numerical instabilities. 
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Therefore, it is more convenient to obtain the expression of ^r,i for a given Ni < (chosen to be iV^ = — ln(l + z^), 
where the initial redshift is Zi — 1.76 x 10''), which gives flr{N = 0) = fJ^.o after solving the differential equation. 

For given £7^,0 and the other parameters we can use Eq. (|27|) to obtain the desired value of ftr^i- At N = Ni, 
Eq. ((77|) is written as 

T , ^K,0 ^2Ni , ^™,0 ^N,\ r> 1 _ 1 ~ ^'".0 ~ - ^r,0 „4(l + s)jV, pl + s foT\ 

i + -p; ^ ~F, ^ I "r,i — ]Tp^ e si^ j . [01) 

Since f2r,i ~ 1 during the radiation domination, it is possible to solve this equation iteratively by assuming that the 
l.h.s. is a small correction (indeed the r.h.s. corresponds to — fiDEi)- 0-th order the solution of Eq. ([3T|l is given 
by 

o(o) ^ I f32l 

1 + (f7K,o/r!,,o) e2^' + (f7™,o/f^r,o) ■ ^ ' 

At first order we find 

(1) _ 1 - (1 - n,n,o - nK,o - O.,o)e4(i+-)^' [l]|."^/r!.,o]i+^ 

1 + (nKfl/nr.o) e2^. + (f7™,o/r!r,o) e^^ • ^ ^ 

This process can be iterated up to the desired precision. In the numerical code, we employ the solution derived after 

(3) 

the three iterations, that is flr^i ~ ^rl- Since the late-time de Sitter background is an attractor, small differences in 
the initial conditions do not lead to very different final solutions. Therefore we indeed find that the three iterations 
are sufficient to derive the parameter Qr,o accurately. 

In the following we first discuss the method for carrying out the likelihood analysis to confront the tracker solution 
with observations and then proceed to constrain the model parameters. 



A. CMB shift parameters 

We use the data of the CMB shift parameters provided by WMAP7, which are related with the positions of the 
CMB acoustic peaks. These quantities are affected by the cosmic expansion history from the decoupling epoch to 
today. The redshift at the decoupling is known by means of the fitting formula of Hu and Sugiyama [40| 

= 1048 [1 -I- 0.00124(f7fc,o/i^)-*'-^3**] [1 + .gi (l^^^o/i^)^^] , (34) 

where gi = 0.0783 {nb^oh^)-°-^^^ /[1+S9.5 (17b,o/i^)°-^^^], 52 = 0.560/[l+21.1 {nb,oh'^)^-^^], h = i/o/[100 kmsec'^ Mpc'^], 
and fib.o corresponds to the today's density parameter of baryons. The shift of the CMB acoustic peaks can be quan- 
tified by the two shift parameters ^l| 



y i2A',o V ^0 H{z)/HqJ Tsiz^) 

where rs{z^) corresponds to the sound horizon at the decoupling, given by 

( \- n 

""^ A. H{z) V3{1 + Silt.o/m^.oi^ + ^)1T ■ 

Note that fl^^ is the today's value of photon energy density and c?i'^'(z*) is the comoving angular distance to the last 
scattering surface defined by d''a\zif) = Tl/[HoyJiT^]. 

The likehhood values of la, TZ, provided by the WMAP7 data [3] are la = 302.09 ± 0.76, TZ = 1.725 ± 0.018, and 
z^ = 1091.3 ± 0.91. The chi-square associated with this measurement is 

la - 302.09 

XcMB = (^a- 302.09, 7^- 1.725, z*- 1091.3) Cci^B I 7^- 1.725 ), (37) 



z* - 1091.3 



where the inverse covariance matrix is 



2.305 29.698 -1.333 
CcMB = I 29.698 6825.27 -113.18 | . (38) 
-1.333 -113.18 3.414 
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B. BAO 

We also employ the data of BAO measured by the SDSS7 The redshift Zd at which the baryons are released 
from the Compton drag of photons is given by the fitting formula of Eisenstein and Hu (42| : 

1291 (fim^o^^)"'^^^ n I m >,2\&2i /qn^ 
= 1 + 0.659 (n^^oh^r^-^ + ) ^ ' (''^ 

where h = 0.313 (r2™,o/i^)"° ''^^[l + 0.607 (a«,o^^)°-^^'^] and &2 = 0.238(rj,„,o/i^)°-^^^- We define the effective BAO 
distance 

Dv{z) = [d\{z) (1 + z)2z/F(z)] , (40) 
where is the diameter distance given by 

dz 



dA{z) = ■== sinh 



H{z)/Ho 



(41) 



The BAO observations constrain the ratio rs{zd)/ Dy{z) at particular redshifts z, where rs{zd) is the sound horizon 
for z — Zd- At z = 0.2 and z = 0.35 the recent observational bounds are rs{zd)/ Dv{0.2) = 0.1905 ± 0.0061 and 
fsizd) / Dv{0-'i5) — 0.1097± 0.0036. The chi-square associated with the BAO is evaluated as 

xIao = irs{zd)/Dv{0.2) - 0.1905, r,(z,)/I?y (0.35) - 0.1097) C^l^ r^ZdyDviO^?^^^^^ ) ' ^^^^ 



where the inverse covariance matrix is 111 



1 -( 30124 -17227 \ 
*-BAo - I -17227 86977 J ' 



C. SN la 



Finally we consider the experimental bounds coming from the observations of the SN la standard candles. The 
apparent magnitudes, together with their absolute magnitudes, can be used to generate the following chi-square (43l | 

2 _ \ " Mobs(Zt) - Mth(^0 

where cr^ ^ are the errors on the data, and /ith is the theoretical distance modulus defined as 

fJ-th{zi) = 51ogio[dL(2:i)] + Mo ■ (45) 
Here dL{z) and are given, respectively, by 

dL{z)^{l + zf HodA{z), Mo = 42.38- 51ogio/i. (46) 
In the following we will make use of the Constitution SN la data sets provided in Ref. [s^ (see also Ref. 13 )• 



D. Observational constraints on the tracker 



We now define the total chi-square as 

= XCMB + xIaO + XSN la • (47) 

In the following we reduce the numerical complexity by fixing several parameters, as ft, = 0.71, q — 0.02258/i^, 
il-y.o = 2.469 X 10~^/i~^, and O^.o = ^7,0(1 + 0.2271 A/'cs) [3], where the relativistic degrees of freedom are set to be 
A/'cff = 3.04. Then two analysis will be performed for the tracker solution: 1) the flat case, VIk,o = 0, for which two 
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0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

s 

Figure 2. Observational contours for the tracker solution in the {s,flm,o) plane for the flat Universe {flK,o = 0). la and 2a 
contours correspond to the internal and external hues, respectively. 



free parameters, flm,o and s, are left to be varied; 2) the non-flat case, for which the additional free parameter, flK,Oy 
is varied as well. 

Later on, when we compare models with different number of free parameters, we will also make use of the Akaike 
Information Criterion (AIC) method (see e.g., [i^). For each model the AIC is defined as 

AIC = xL„ + 27', (48) 

where V is the number of free parameters in the model, and Xmin minimum value for in the chosen parameter 

space. The smaller the AIC the better the model. To be more precise, if the difference of between two different 
models is in the range < A(AIC) < 2, the models are considered to be equivalent, whereas if A(AIC) > 2, the data 
prefer one model with respect to the other. 



1. Flat case: ^1k,o = 

In this case we compute the on a grid in the intervals < s < 0.9 and 0.25 < flrn,o < 0.32. The minimum value 
of is found to be Xmin ~ 468.876 for the model parameters s — 0.03446 and il.m,Q — 0.27159. Then we calculate the 
difference of x^ at each grid point, that is, Ax^ = ~ Xmin- When Ax^ > 2.88 the chi-square distribution, with two 
free parameters, excludes the models with those values of x^ at 68% confidence level (Icr), whereas when Ax^ > 5.99 
those models are excluded at 95% CL {2a). 

Our numerical results are plotted in Fig. [21 Even if we use the Gaussian likelihood function P cx e^^^^^, we find 
that the observational contours are similar to those given in Fig. [21 The parameters s and flm,o are constrained to be 

s = 0.034t^;327, a„^o - 0.271™;^ (95% CL). (49) 

This shows that the tracker solution with —1.36 < wde < — 1 during the matter era can be allowed observationally. 
The ACDM model corresponds to the line s = in Fig. [2l which, as expected, is inside the la contour for 0.264 < 
^m,o < 0.273 (with Xacdm ~ 469.024). The best-fit x^ for the extended Galileon model is slightly smaller than that 
in the ACDM model. However, since A(AIC) = 1.85, the observational data do not particularly favor the extended 
Galileon model over the ACDM model. 
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Figure 3. Observational contours for the tracker solution in the {s,Qm,o,^K,o) space with the cosmic curvature taken into 
account, la and 2a contours correspond to the black and red regions, respectively. 



2. Non-flat case: Qk.o 

In the presence of the cosmic curvature K we also evaluate the in the parameter space given by < s < 0.9, 
0.25 < flm,o < 0.32, and -0.006 < ilxfi < 0.014. The minimum value of is found on the point s = 0.067, 
rimfi = 0.2768, and flxfi = 0.00768, at which xLin,K = 467.436. We evaluate the difference Axk = - Xmin,K 
on each grid point, according to which we can exclude: a) the values of at 68 % CL (which, for a system with 3 
parameters, corresponds to > 3.51), and b) the values of at 95% CL (i.e. ^Xk > 7.82). In Fig. [3] we plot a 

three dimensional region (by showing only some slices of it) of the allowed parameter space for the tracker solution. 
We find that the model parameters are constrained to be 

s = 0.067tS];3|3 ^ l]™,o = 0.277t°;°22, rJ/^.o = 0.0077^^9 (95% CL). (50) 

The AIC for the best-fit extended Galileon model is AIC = 473.436, whereas in the best-fit ACDM model with 
two parameters ilm,o and fi_R-,o the AIC is found to be AICacdm,^ = 472.543 with Xacdm k ~ 468.543. Since 
A(AIC) = 0.893 between the two models, they are equivalently supported by the data. However, if we compare the 
non-flat extended Galileon model with the flat ACDM (1 parameter only, 51™, o), then we see that A(AIC) = 2.412, 
which states that the flat ACDM is better supported by the data. According to this criterion, the flat extended 
Galileon tracker is favored over the non-flat case from a statistical point of view. 



IV. COSMOLOGICAL PERTURBATIONS AND THE DISCRIMINATION BETWEEN THE MODELS 

The property of the tracker solution does not depend on the parameters a and /3 at the level of the background 
cosmology. If we consider cosmological perturbations, it is possible to discriminate between the models with different 
values of a and /3. In order to confront the extended Galileon models with the observations of LSS, CMB, and weak 
lensing, we shall study the evolution of matter density perturbations as well as gravitational potentials. Since our 
interest is the growth of non-relativistic matter perturbations in the late Universe, we do not take into account the 
radiation as far as the perturbations are concerned. 
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A. Linear perturbation equations 



Let us consider the perturbed metric in the longitudinal gauge about the flat FLRW background [4£ 

ds^ = -(1 + 2*) dt^ + a^{t){l + 2$)(5yda;*dx^ , 



(51) 



where ^ and $ are scalar metric perturbations. We perturb the scalar field as 4){t) + 6(l){t,x), and non-relativistic 
matter as well, in terms of the matter density perturbation Sp„i and the scalar part of the fluid velocity v. The density 
contrast of non-relativistic matter is defined as ^ = Spm/ Pm- We also introduce the gauge-invariant density contrast 



(5„. = 6 + 3Hv . 



(52) 



The full linear perturbation equations in the Horndeski's most general scalar-tensor theories were derived in Ref. [47 
For the extended Galileon models, in Fourier space, they are given by 



/j2 p 



- p„i(5 = , 



Be* 

^3$ 



AeS(l)- Ai-^/i 

3^6$ + D2S(j) + 4 

(5 + 3$+— w = 0. 



• + PrnV = , 



Ao^ 



^2 



A. 



fc2 



^11 * = 0: 



where k is a comoving wave number, and 

Ai = 12HG4 - 6 ^XG3,x - 48 HX (G4,x + XG^^xx) - QH^X (5 G^.x + 2 XG^.xx) <i> 
- - iK,x + 2 XK.xx) 4> - (^HX (3G3,x + ^XG^.xx) - 6H^' (3 Gi,x + 12 XG^^xx 
- 2H^X (15 Gs.x + 20XG5,xx + 4X2G5,xxx) , 
A3 = 4G4 - SXGi^x - 4:H(pXG5,x , 

2X (X X + 2XK,xx) - UH^Gi + p„, + 12i/X (2G3,x + XG3,xx) ^ 



AX^G^ 



A.XXX I 



A4 



+UH^X (7G4^x + l6XG4,xx + ^X^G^^xxx) + ^H^X {ihG^^x + HXG^^xx + 2X^G^^xxx) 
As = -2XG3.X - 4.H (G4^x + 2XGi^xx) 4> - 2H^X {SG^^x + 2XG5,xx) , 
Be = 4G4 - 4XG, 



5,X(i 



Br = -AGi.xH4> - 4(G4,x + 2XGi,xx)4> - ^{G^,x 



H). 



(53) 
(54) 
(55) 

(56) 

(57) 

(58) 

(59) 

(60) 
(61) 

(62) 
(63) 
(64) 
(65) 
(66) 



2H) 

XG5,xx)^. (67) 



XG5,xx)Hcj)cb~'iXG5,xiH 
K,x^ + QHXG^^x + &H\Gi,x + 2XGi,xx)<P + 2XH''{iGz,x + 2XG5,xx) , 
—K x — 2(G3^x + XG3^xx)4' ~ '^HG3,x(t> — AHiZG/^^xx + 2XGa 

,xxx)4'4' ~ 2G4,x{'iH 

-AXGi,xx{5H^ + 2H) - 2H^G5,x + ^XG^^xx + 2X^G5,xxx)4> - '^HiH^ + H){G^,x 

The readers may refer to the paper [47] for the explicit forms of the coefficients D2, D3, D4, Ds, and Du. Since the 
perturbation equations (|53|) -(|58 p are not independent, we do not need to know these unwritten coefficients to solve 
the equations numerically. Now we are dealing with a massless scalar field, so that the mass term M does not appear 
in the perturbation equations. 

From Eqs. (|57|) and ([55)1 the gauge- invariant matter perturbation ([5^ obeys 

fc2 



* = 3 / + 2HI 



S,n + 2HS„ 

where / = Hv — We define the effective gravitational potential 

$eff = («- - *)/2 . 



(68) 



(69) 



This quantity is related to the deviation of the light rays in CMB and weak lensing observations [48j. To quantify the 
difference between the two gravitational potentials <i> and 5*, we also introduce 



by which Eq. can be written as $off = *(1 + v)/2- 



(70) 
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B. Quasi-static approximation on sub- horizon scales 



For the modes deep inside the Hubble radius (}? ja^ ^ i/^) we can employ the quasi-static approximation under 
which the dominant contributions in the perturbation equations are those including j and (5 (or (49l.[50|. This 
approximation is known to be trustable as long as the oscillating mode of the field perturbation is negligible relative 
to the matter- induced mode. Combining Eqs. (|53l) . ([M)) . and ([55]) under this approximation, it follows that [l^l 

^ ~ -47rGoff/9™(5 , (71) 



a2 



where 



Here we introduced the bare gravitational constant G = l/(87rMpj). Substituting Eq. (|7ip into Eq. (I68p with the 
relation S„i ~ 5 (valid for k"^ /a^ ^ -ff^), we obtain 

C + (2 + C - l^^raSra ^ , (73) 

where a prime represents a derivative with respect to TV = In a. 

Under the quasi-static approximation the quantity rj defined in Eq. (|70p reads (47] | 



On using Eq. ((7T|) . the effective gravitational potential <i>cff = ^(1 + '/)/2 yields 



3, f aH^ ^ 



^^s'^-^iy—j f^mf^m, (75) 
where 

^~ G 2 ~ AlBe + AlDg - 2A3AeB7 ' ^ ' 

The ACDM model corresponds to K = -A, G4 = A/pi/2, 63^65== 0, in which case A3 = = 2M'^^, = 0, and 
By = 0. Then one has G^s/G = 1 and ^ = 1 from Eqs. ([7^ and (jTS]) . 

In the extended Galileon models the general expressions of G^s I G and ^ are quite complicated. In what follows we 
shall focus on the evolution of cosmological perturbations for the tracker solution (r\ = 1). 

In the early cosmological epoch (r-i <C 1) and during the matter domination (51,. = 0), Gcn/G and ^ are approxi- 
mately given by 

Geff/G~ 1 + [27p(2p- l)(3a2p + 6/32(2p- 1) + /?(1 + 3a - 2(1 + 6a)p)) + 9(6aVl8p - 5) 

+2/3(2p - l)(/3(90p - 3) - lip) + a{2p{2p - 1) - 3/3{l + 4p(33p - I3))))q + 2(-9(9a - 23/3){a - 2(3) 
+&{U7a^ + a(5 - bQip) + 6/3(71/? - 3))p -f 8(9a - 18/? - 1)^^ + 16p3)(?2 + 4(9(27a - 46/3)(a - 2^) 
-12(9a - 18;9 - 8)(a - 2/3)p + 16(1 - 3a + QI3)p^)q^ + 48(a - 2/3)(6/3(7 - 2p) - 8p + 3a(2p - 7))q^ 
+576(a-2/3)\>2/A, (77) 
e ~ 1-1- [27p(2p - l){Za^p + 6/3^(2p - 1) + ^(1 + 3a - 2(1 + 6a)p)) -I- 9(24a^p(4p - 1) 

+2p{2p - l)(/3(78p 3) - lip) + a(4p(2p - 1) + /3(3 + 12(10 - 29p)p)))q + 2(-45(a - 2^)/3 
+6(96a^ + a(10 - 321/3) -I- 6/3(43/3 - 2))p + 4(27a - 54/3 - 4)p2 -f 32p3)q2 ^ 8(9(30, _ io/3)(a - 2/3) 
-3(27a - 54^ - 26)(a - 2f3)p + 16(1 - 3a + 6P)p'^)q^ + 48(a - 2/3)(66^ - 8(2 + 3l3)p + 3a(4p - ll))g'* 
+1152(a-2/3)2q5]r2/(2A), (78) 

where 

A = 4g[(3 - 9a)p + 12p^ + 2pq{a{60 - 66q) + 20q - 11) + 2p^{a{9 - I8q) + 22q - 6) - 3ag(9 - 54g + 40g^) 

+3/3(2p + Aq- 1)(3 - 22q + 20q'^ + Qp{2q - 1))] . (79) 
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Figure 4. Allowed parameter space in the (a,/3) plane for p = \ and q — 5/2. In the blue region the conditions for the 
avoidance of ghosts and Laplacian instabihties of scalar and tensor perturbations are satisfied. The right panel shows the 
enlarged version of the left panel in the (a, /3 — a/2) plane. The borders correspond to (a) /3 = a/2, (b) /3 = a/2 — 1/15, (c) 
/? = (408q + 68 - 2\/T7v/3(272 - 75a)a + 68)/561, and (d) /3 = (242 - 15a + 4V3630 - 495a)/99, respectively. Taken from 
Ref. [iil. 



To derive Eqs. ([77)1 and ([75)1 we set 57,, — and performed the Taylor expansion around r2 — 0. 
At the de Sitter solution (ri = ?'2 = !)■ and ^ are simply given by 

Gcff/G = e 



2(l-j5)+3(l + 2g)(a-2/3) 



(80) 



The equality of Goff / G and ^ comes from the fact that 77 = 1 at ri = r2 = 1 . 

Let us first consider the theory where a = /3 = 0. From Eqs. ([77|) - ([50| it follows that 



Geff/G = e - 1 

Gcff/G = ^ = 



Apq 



6p + IO9 - 3 



r2 



1 



1-p 



(for r2 < 1), 
(for r2 = 1). 



(81) 
(82) 



For p — 1 [30] both Gcff/G and ^ diverge at the de Sitter solution. In this case one has Gcff/G = £. — l + 4:qr2/ {lOq + 3) 
in the regime r2 <C 1, so that Gcff/G and ^ are larger than 1 for g > 0. The property that ^ is as large as Gcff/G 
leads to the enhancement of the effective gravitational potential relative to the matter perturbation dm normalized 
by a |;30]. This gives rise to the anti-correlation between the late-time ISW effect and the LSS. Then the parameter q 
is constrained to he q > 4.2 x lO'^ at the 95 % confidence level |33], which means that wde is very close to —1 along 
the tracker. 

The situation is different for a 7^ and /3 7^ 0. Let us consider the models with p — 1 and q = 5/2, i.e. s = 0.2, in 
which case the models are compatible with the observational constraints discussed in Sec. IIIIl In the regime r2 ^ 1, 
Eqs. dZZl) and ^ read 



1 



Gcff/G 

e^l 



r2 , 



48411^2 _ 3a(3560 + 60291/3) + 22(50 + 924/3 + 7641/32) 

10(308- 1536a + 3201/3) 
75276a2 _ 3a(8020 + 99981/?) + 22[100 -I- 3/3(733 -t- 4527/3)] 



20(308 - 1536a -I- 3201/3) 



-r2 ' 



At the de Sitter solution Eq. ([50]) gives 



Geff/G — ^ 



1 



9{a - 2/3) 



(83) 
(84) 

(85) 



In Ref. [36| the authors clarified the viable parameter region in which the ghosts and Laplacian instabilities of 
scalar and tensor perturbations are absent. In Fig. |3] we plot the allowed parameter space in the {a, f3) plane for 
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Figure 5. The left and right panels show the values of Gch/G and ^ = (Gcff/G)(l + t?)/2 versus a for p = 1 and q = 5/2, 
respectively, along the tracker solution (ri = 1) at r2 = 0.1. The blue regions illustrate the viable parameter spaces in which 
the parameters a and f} belong to the blue region in Fig. [l] The borders (a), (b), (c), (d) correspond to those given in Fig. [l] 



The viable region is surrounded by the four borders: (a) [3 — a/2, (b) j3 
17^3(272 



p = 1 and (7 = 5/2. 

P = (408a + 68 - 2%/T7V3(272 - 75a)a + 68)/561, and (d) ^ = (242 - 15a + 4^/363(^ 
a — P — Q \s on the border lines [the intersection of the lines (a) and (c)]. in which case both Gcs/G and ^ diverge 
at the de Sitter solution. 



a/2 - 1/15, (c) 
495^)/99. The model with 



Figure [S] illustrates the regions for the possible values of Gcs/G in Eq. ([55)1 and ^ in Eq. for p = 1 and q — b/2 
at r2 — 0.1. The blue regions are surrounded by the borders (a), (b), (c), (d), which correspond to those given in 
Fig. 0] respectively. The lines (a) in Fig. [S] show Geff/G and ^ for /? = a/2, along which both Gcs/G and ^ decrease 
for larger a. One has Ges/G = ^ ~ 1.035 for a = /3 = 0, in which case our numerical simulations show that there is 
an anti-correlation between Sm/a and $cff for the modes deep inside the horizon. When a > the variable ^ is larger 
than Geff/G on the line (a) for the same a, so that the anti-correlation tends to be even stronger than that for a = 0. 

On the line (b) in Fig.[4l i.e. (3 — a/2 — 1/15, one has Ges/G ~ 1 -|-0.52r2 and ^ = l-|-0.16r2, which are independent 
of the values of a. The fact that Gcs/G is larger than ^ may imply the absence of the anti-correlation between 5m /a 
and $cff- In fact we will show in Sec. IIV CI that the anti-correlation tends to disappear as the model parameters 
approach the border line (b) in Fig. S) For given a, as we move from /3 = a/2 to /3 = a/2 — 1/15, G^s/G increases 
whereas ^ gets smaller. For < a < 34/75 the viable regions for Gcs/G and ^ in Fig. [5] are surrounded by the lines 
(a), (c). and a = 34/75, whereas for 5.579 < a < 5.646 they are surrounded by the lines (b), (d), and a = 5.579. 

In Fig.[6]we plot G^s/G = £, = l/[9(a — 2/3)] for the late-time de Sitter solution. Both Goff/G and ^ diverge on the 
line (a) in Fig. [H including the case a — (3 — Q. On the line (b) one has Gcs/G = ^ = 5/6, so that the growth rate 
of matter perturbations is smaller than that in the ACDM model around the future de Sitter solution. The allowed 
values of Gcs/G and ^ exist in the wide regions in Fig. [51 

The present epoch corresponds to the regime between r2 ^ 1 and r2 = 1, so we need to resort to numerical 
simulations to estimate the growth rate of perturbations accurately. 



C. Numerical simulations 



In order to study the evolution of perturbations for a number of different wave numbers (both sub-horizon and 
super-horizon modes), we shall solve the full perturbation equations numerically without using the quasi-static ap- 
proximation. Let us introduce the following dimensionless variables 



V = Hv, S(p = 6(l)/{xdsMpi) , 



(86) 
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Figure 6. Ges/G versus a for p = 1 and g = 5/2 at the de Sitter solution (ri = 1, r2 = 
the same as Gbh/G. Goff/G diverges at /3 = a/2, so that it is unbounded from above, 
those given in Fig. [l] respectively. 



1). Note that the parameter ^ is exactly 
The borders (b), (c), (d) correspond to 



with ii = A,/{HM^^), A2 = Xds^2/(i?Mpi), ^3 = ^a/Mp^i, A4 = AJ{HH-II,), A^ = XdsA/Mpi, = B^/M^^, 
Br = XdsBj/Mpi, and C4 = xasCa/ [HM^^). From Eqs. and ^ we obtain 

^ ^-{Be<i> + Bj6ip)/Az, (87) 
= [(314^6^6 + iii2^6 - 3i2igA:V(ai^)' - ^A^B^nra)<^ 

+ (312^3^4 + iAiAeB^ + iiiz-By - gie^yO^ - ZA^Alk'^ / {aHf)5v + giaig^^miJ + 9^2^^^™^] 

X [313(^1^6-^2^3)]"', (88) 

Sif' = -[(i?i?6 + 3^3^-66 - iAlk^l{aHf - g^iJe^^m)^' 

+ {AlBr + 3i3^i?7 + 3iii3C'4 - 3ip6fcV(ai?)^ - 9i3i?7f^™)'5(y5 + 9i^f7,„5 + 9iii3f^m^] 
X[3i3(ili6-i2i3)]"', (89) 

5' = -W -e/{aHfV , (90) 
= (i77i/)y + ^' . (91) 

In order to recover the General Relativistic behavior in the early cosmological epoch we choose the initial conditions 
$' = 0, 5if' = 0, (5iy9 = 0, and 5 — 10^^, in which case $i, ^f^, and Vi are known from Eq. (|87p - (|5^ (where the subscript 
"i" represents the initial values). For non-zero initial values of 5lp the field perturbation oscillates at the early stage. 
For the initial conditions with \5(fi\ the evolution of perturbations in the low-redshift regime is hardly affected 

by the oscillations. This situation is similar to that found in Ref. [s^l for the model with a — [3 — 0. 

In Fig. [7| we plot the evolution of (5m /a and <i>off versus the redshift z — 1/a— 1 forp = 1, g = 5/2, a = 3, and 
(3 = 1.49 with several different wave numbers. In this case the model parameters are close to the border line (a) in 
Fig. m The soHd curve in Fig. [7] corresponds to the simulation for the wave number k — 300aoi?o — 0.1 hMpc~^, 
where the subscript "0" represents the today's values. The numerical results of Sm/a and show good agreement 
with those derived under the quasi-static approximation on sub-horizon scales. For the model parameters used in 
Fig.[7|the analytic estimates and give Gcs/G ~ 1 — 0.05r2 and ^ ~ 1-1-0. 23r2 in the regime r2 <C 1, whereas 
from Eq. ([55)1 one has Gcff/G — ^ — 5.56 at the de Sitter solution. From Fig. [7|we find that Sm/a is anti-correlated 
with $cff for the modes k ^ floiJo- This property comes from the fact that the parameter ^ is always larger than 1, 
while Goff/G is smaller than ^ in the regime r2 ^ 1. In Fig.[7]the growth rates of 6m and $off decrease for smaller k. 
The anti-correlation between Sm/a and $off is present for the modes k > 5aoi?o- 

Figure [5] shows the evolution of 6m/o, and $eff for p = 1, q = 5/2, a = 3, and /? = 1.45 with several different wave 
numbers. In this case /3 is smaller than that used in Fig. [T) From Eqs. (|83|) and (|84|) one has Geff/G = 1-1- 0.25r2 
and ^ = 1 + 0.16r2 for r2 <C 1, whereas Gcg — £, — 1-11 at the de Sitter solution. Compared to the case (3 = 1.49, 
Gcff/G and ^ get larger and smaller, respectively, for the same value of r2 (^ 1). This leads to the suppression of the 
growth of $cff- From Fig. [5] we find that Sm/a and $cff are positively correlated for the mode k = 300ao-ffo- For the 
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Figure 7. Evolution of 5m/a (left) and $cff (right) versus the redshift z for p = 1, g = 5/2, a = 3, /3 = 1.49 along the tracker 
solution. Note that hraja and <l?cft are normalized by their initial values, respectively. Each curve corresponds to the evolution 
of perturbations for the wave numbers k — 300ao-?/o, k — 15aoHo, k — 5aoHo, and k — aoHo- The dotted curves show the 
results obtained under the quasi-static approximation on sub-horizon scales. In this case there is an anti-correlation between 
Sm/a and $efi for the modes k > 5aoHo. 




Figure 8. Similar to Fig. [T] but for the model parameters p = 1, q — 5/2, a = 3, /? = 1.45. In this case 5m/a and <E>cfi are not 
anti-correlated. 



wave numbers k ^ aoHo the quasi-static approximation reproduces the numerical results in high accuracy. On larger 
scales both 5m and $off evolve more slowly, so that Sm/a and $eff are also positively correlated. 
We define growth index 7 of matter perturbations, as [5l| 

^^{nmr- (92) 

In the ACDM model 7 ~ 0.55 for the redshift < z < 1 [52|. In Fig. [U] we plot the evolution of 7 for p = 1 and 
q = 5/2 with the wave number k = 300aoi?o in three different cases: (i) a = 3, /? = 1.49, (ii) a = 3, /3 = 1.45, and 
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Figure 9. Evolution of the growth index 7 for the mode k — SOOaoHo in the regime < 2: < 1 with the model parameters (i) 
p=l,q = 5/2, a = 3, P = 1.49, (ii) p=l, q = 5/2, a = 3, /3 = 1.45, and (iii) p = 1, g = 5/2, a = 3, /3 = 1.434. 



(iii) a = 3, /? = 1.434. The numerical simulations for (i) and (ii) correspond to the model parameters used in Figs. [7] 
and |51 respectively, whereas the model parameter in the case (iii) is close to the border (b) in Fig. [H In the cases 
(i), (ii), (iii) the effective gravitational couplings in the regime r2 ^ 1 are Gcs/G = 1 — 0.05r2, Gcs/G = 1 + 0.26r2, 
Ges/G = 1 + 0.51r2, respectively, whereas at the de Sitter solution Gcs/G = 5.56, Ges/G — 1.11, Gcs/G — 0.84, 
respectively. For smaller /3, GcS is larger in the early epoch (r2 <C 1), so that the deviation of 7 from the value 0.55 
is more significant. On the other hand, Gcff at the de Sitter solution gets smaller for decreasing /3, which leads to the 
approach to the value 0.55 around z — 0. 

While the above discussion corresponds to the case s = 0.2, we have also studied the evolution of perturbations 
for different values of s and found the similar properties to those discussed above. For p = 1 the upper bound for 
the allowed parameter space in the {a, /3) plane is characterized by the line /3 — a/2 [s^, around which the ISW-LSS 
anti-correlation is present for < s < 0.36 (i.e. for s constrained observationally at the background level). As /3 gets 
smaller for a fixed a (> 0), the anti-correlation tends to disappear. We also found that, for p = 1 and < s < 0.36, 
the growth index 7 varies in the range 0.4 < 7 < 0.6 for < z < 1. The ISW-LSS anti-correlation as well the variation 
of 7 allows us to discriminate between the extended Galileon models and the ACDM. 



V. CONCLUSIONS 



In this paper we have studied cosmological constraints on the extended Galileon models of dark energy. For the 
functions ([7]) with the powers there exist tracker solutions along which the field equation of state wde changes as 
— 1 — 4s/3 (radiation era) — ^ — 1 — s (matter era) — ^ — 1 (de Sitter era), where s = p/ (2g) > 0. Unlike the case of the 
covariant Galileon (s = 1), wde can be close to —1 during the radiation and matter eras for < s ^ 1. Moreover, 
even for wde < ~1, there are viable model parameter spaces in which the ghosts and Laplacian instabilities are 
absent. 

Using the recent data of SN la, CMB, and BAO, we placed observational constraints on the background tracker 
solutions with s > 0. In the flat Universe we found that the model parameters are constrained to be s = 0.034lg'Q34 
and flm.o — 0.2711q'o^q (95% CL) from the joint data analysis. The chi-square for the best-fit case (s = 0.034) 
is slightly smaller than that in ACDM. However the difference of the AIC information criteria between the two 
models is A(AIC) = 1.85, so that the extended Galileon is not particularly favored over the AGDM. We also carried 
out the likelihood analysis in the presence of the cosmic curvature K and obtained the bounds s = 0.067to o67j 
n„„Q = 0.277l^;J522, and f^^.o = O.OOTJtoMlr (95% GL). The tracker for the covariant Galileon (s = 1) is disfavored 
from the data, in which case only the late-time tracking solution is allowed observationally [3ll|. 

The background quantities for the tracker are independent of the values of a and (3. This means that the models 
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with different a and /? cannot be distinguished from the observational constraints derived from the background cosmic 
expansion history. In order to break this degeneracy we studied the evolution of cosmological perturbations in the 
presence of non-relativistic matter for the flat FLRW background. As shown in Ref. [3d|, for a = /? = 0, the matter 
density perturbation divided by the scale factor a is anti-correlated with the effective gravitational potential $off 
for the modes relevant to the LSS. This leads to the anti-correlation between the LSS and the ISW effect in CMB, so 
that the parameter s for the tracker is severely constrained to be s < 1.2 x 10^'' (95% CL). 

For the models with a ^ and /3 7^ 0, however, the correlation between 5m/ o- and <i>cff depends on the values of a 
and /3. If p = 1 and q = 5/2 (i.e. s = 0.2), for example, Sm/a and $cff tend to be positively correlated for the model 
parameters close to the border (b) in Fig. |l](/3 — a/2 — 1/15), whereas they show anti-correlations for a and /? close 
to the border (a) (/? = a/2). The typical examples of the positive and negative correlations are plotted in Figs. [5] 
and [71 respectively. The qualitative differences between these two cases can be understood by estimating the effective 
gravitational coupling Goff and the quantity ^ — (Gcff/G)(l -I- 77)/2 derived under the quasi-static approximation on 
sub-horizon scales. As the model parameters approach the border (a) in Fig. [4l ^ gets larger while GoS decreases, so 
that the anti-correlation between 5m/a and 3>eff tends to be stronger. We studied the evolution of perturbations for 
different values of s in the range < s < 0.36 and found that the basic properties for the ISW-LSS correlation are 
similar to those discussed for s = 0.2. 

We also estimated the growth index 7 of the matter perturbation and found that, for p — 1 and < s < 0.36, it 
typically varies in the range 0.4 < 7 < 0.6 at the redshifts for < z < 1. Hence it is also possible to distinguish 
between the extended Galileon models and the ACDM model from the galaxy clustering. However, we expect that 
the tightest observational bounds on the values a and (3 should come from the ISW-LSS correlation. We leave such 
observational constraints for future works. 
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